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DESIGN PFOCEDURE FOR SATISFYING TIME DOMAIN BOUNDS 
FOR NONMINIMUM-PHASE SYSTEMS 

Abstract— This paper presents two design techniques 
applicable to nonminimum-phase systems. Both are de- 
signed to handle plants with one right-half-plane zero 
which may vary and any other variation of the plant 
parameters within known limits. The specifications 
that must be designed to are given as a set of step 
response bounds in the time domain. A completed de- 
sign will yield responses that stay within the time 
domain bounds at all times and utilize the entire 
region of allowed variation. 

The first method is an analytic one which finds an 
equivalent set of minimum-phase step response bounds. 
Once these are found, ciirrent design techniques may be 
used to complete the design. This method can be ex- 
tended easily to design for multiple right-half-plane 
zeros, but the actual design will be very complicated 
for this probl^. When there is variation of the 
right-half^plane zero, this method requires many de- 
signs before the best one is found. For these reasons, 
this method is not very useful. 



The second method makes a direct transfer of the 
time domain hounds to the frequency domain. From these 
new boundaries and other resulting specifications, a 
worst case zero is found that will yield a design that 
is satisfactory as the zero varies. This method is 
thoroughly derived and explained for a plant with one 
right-half-plane zero and the plant does not depend on 
the position of the zero in any other way. These de- 
sign techniques are the first contributions to the non- 
minimum-phase problem. 
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CHAPTER I 


INTRODUCTION 
1.1 Problem Statement 

The purpose of this paper is to present two methods 
of design which are applicable to nonminimum-phase prob- 
lems. Both methods are thoroughly derived for one right- 
h®lf~pl3J^® zero and the extensions are indicated in the 
first method to handle multiple right— half— plane zeros. 
One method of design is analytical whereas the other is 
best suited to digital computer implementation. Both 
methods do permit insight into the means of design. The 
program utilized in the computer implementation of the 
second example is in the appendix. 

The problem specifications are given in general as 
a set of time domain boundaries within which the final 
step response must remain at all times. There may be 
ignorance in the right— half— plane zero(s) and in any 
other parameters the plant may contain. The only thing 
known about the parameters are the limits that they are 
permitted to vary within. The rate of pareimeter varia- 
tion is assumed to be slower than the system response 
time so time dependence of the parameters may be neg- 
lected. A successful design will yield satisfactory step 
responses for all possible sets of plant parameters. 



The step response boundaries are unique because 
they restrict the initial undershoot that must occur in 
a nonminimum-phase system. The optimum design will per- 
mit step responses that utilize the entire allowed re- 
gion of variation, both negative and positive. In this 
way the compensation is not more than it needs to be in 
complexity and bandwidth. 

The feedback structure used in the design is a two 
degree of freedom structure shown in figure 1.1. A two 
degree of freedom structure is used because both sensi- 
tivity and minimum bandwidth must be designed for. The 
feedback structure is required to meet the sensitivity 
constraints, i.e., the plant variation cannot result in 
more variation in the step response then allowed by the 
step response boundaries. The prefilter F will adjust 
the bandwidth to the minimum allowed so as to make the 
system as immune to noise as possible. 
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Figure 1.1 CloBed Loop Structure 
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1 . 2 History of the Problem 

Right-half -plane zeros occur in many circumstances. 
Some examples are where there are multiple channels over 
which information can be passed, lattice networks, trans- 
mission lines, piping-flow delays, transportation lags, 
semiconductor diffusion tU # heat exchangers, mercury 
thermometers [2] , and in aircraft control [3] . The pre- 
sence of a right-half-plane zero greatly affects the 
stability conditions of a design. The loop transmission 
cannot be allowed to decrease too quickly or the addi- 
tional phase lag will be so great as to cause condition- 
al or complete instability. In a stable nonminimum- 
phase system, there is an upper limit to the crossover 
frequency [4] . A nonminimum-phase system is one that 
has a zero in the right half of the s-plane. Nonminimum- 
phase systems are slow in response because they initial- 
ly go in the wrong direction [5] . A magnitude plot will 
completely determine the transfer function of a minimum- 
phase system, but both magnitude and phase plots are re- 
quired for a nonminiraum-phase system. 

These considerations make the design for a non- 
minimum-phase plant a unique one, and from the areas 
where the problem arises it is an important one. In- 
spite of this, there is no design technique available 
for the solution of this problem. 
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1 . 3 Method of Approach 

There is a design procedure presented in a paper 
written by Horowitz and Sidi which can formulate a de- 
sign for a plant with parameter variations from a set 
of minimum-phase time domain bounds [6] . The first method 
presented here solves for the set of minimum-phase bounds 
equivalent to the given set of nonminimum- phase bounds. 
That is, a design meeting the minimum-phase bounds will 
also meet the nonminimum-phase ones. The minimum-phase 
bounds are found by solving the differential equation 
relating them to the nonminimum-phase bounds. Because the 
goal of this method is to yield an equivalent set of 
minimum-phase time domain specifications, this method is 
called the "time domain method.” The remainder of the 
design is completed by the above procedure. 

The second method approximates the transfer function 
by a third order function and finds a set of boundaries 
on the allowed variation of the magnitude of the transfer 
function in the frequency domain. This is done by taking 
the composite of the Bode plots of many transfer func- 
tions which have satisfactory step responses. It was 
found that a design for the smallest right-half -plane 
zero and the lowest set of specifications found above 
will be satisfactory as the right-half -plane zero 
varies. Since the design is carried out for a fixed po- 
sition of the zero, a magnitude plot will be sufficient 
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to describe the transfer function. Thus, once the trans- 
formation to the frequency domain is carried out, the 
Horowitz and Sidi design technique can be used to com- 
plete the design. 



G 

CHAPTER II 
TIME DOMAIN METHOD 

2 . 1 Motivation of Method 

Horowitz and Sidi have written a paper [6] which 
presents a design method utilizing a set of minimum-> 
phase step response boundaries. One design approach for 
a set of nonminimum-phase time domain boundaries is to 
start by finding an equivalent set of minim\im~phase 
tolerances and applying the Horowitz and Sidi method to 
these tolerances to complete the design. This transfor- 
mation can be made by solving the differential equation 
relating the nonminimum-phase (nmp) transfer function 
and its minimum-phase (mp) component. This method will 
be known as "the time domain method" because its goal is 
to yield an equivalent set of mp time domain specifica- 
tions and is not concerned with the techniques used to 
complete the design once the equivalent bounds are found. 

2.2 Technique 

If there is only one right-half-plane (RHP) zero, 
the transfer function may be written as 

T(S) = T^(s) (1 - |) 

where T^(s) is the minimum-phase portion of T(s), and 
"a" is the location of the RHP zero. If the input is 
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R(s) = — 0 a unit step, then the output response C(s), 
assigning zero initial conditions, can be written as 


C(s) = R(s)T(s) = = -i2 E 

s s 


c{t) = c (t) - ~ 6 (t) 

m am 


where the roles of R(s), C(s), and T(s) are clarified 
in figure 2.1 and 




-1 


m 


This means; 




C(S) - i 


C (0) 
m 


i-f 


( 2 . 1 ) 


for the general case when there are initial conditions. 

This equation points out that an initial condition 
is required to solve for the mp set of step response 
bounds. A point is required where both the mp and nmp 
responses are known or the same so the initial condition 
can be solved for. This is a point where the response 
has settled to a constant. Start at this point and work 
backwards solving for the initial conditions. An ini- 
tial condition is required at every time when the nmp 
boundary is discontinuous. 
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2 . 3 Example of Design Method 

To demonstrate the technique the equivalent rap 
boundaries for the nmp ones in figure 2.1 will be found. 
The results are shown dotted in figure 2.1, the mp 
bounds. From t 2 to t^, the upper bound is constant at 
1 + a. The equivalent mp boundary can be found by solv- 
ing equation 2,1. The initial condition at t 2 can be 
solved for by knowing that at t^ both the mp and nmp 
boundaries must equal one. From equation 2.1; 


m 


l+g _ 1 
s "a 


m 


(0) 


1 


£ 

a 



Figure 2,1 Hypothetical set of nonminimum-phase bound- 


aries and equivalent set of minimum-phase 


bounds . 
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which yields 


C (t) =l+a- (1+a-C (t,) 
lu m ^ 


To solve for , the initial condition, one can 

solve the equation at a point where Cj^(t) is known, t^, 

= l- l + a- (1 + a- Cj^(t2))e^^4- ^ 3)3 


which yields 


C (t„) = 1 + a - ae 
in z 




Thus the value at t 2 is known and can be used to solve 
for the value at t = 0 and thus specify the entire upper 
bound. For this region; 

a 


- 

m 


1 - ^ 
a 


which yields 


C_(t) 

m 


I ^ - <1 - 


C (t_) = 1 + a - ae 
m z 


-a(t4-t2) 


d d 


r "^(t.-t,,) , -at- 

C^(0) = [1 + a - ae * ^ . <| + dt^jje ^ ^ | 

Now the entire equivalent mp upper boundary is known. 

The same method is utilized to find the mp lower 
bound. Starting at t^ as before: 
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1-6 


C^(s) = 
m 




C^(t) ^1-B - (1-6 - Ci^<0))e 


a (t-t^) 


'=m<'^4' -1 = 1-6- (1 - 6 - C^(t3>)e 




C (t,) = 1 - 6 + 6e 
m J 


-aCt^-t^) 


now for < t < 


C(t) = dt - A 

a a <=m"» 


.2 s 


c (s) = 
m 


1 - ^ 
a 


d d a(t-t,) 

C^(t) = I - A + dt + [A - I + C^(t^)]e 


C (t-) = 1 - 6 + Be 
m 3 




d d 

= ^ - A + dt3 + [A - - + C^(t^)]e 


C (t, ) = - - A + 
in 1 a 


A - - + Be ^ ^ 

a 


-a(t3~t^) 


for the boundary from 0 < t < t. 


a C-n'O' 


C (s) = 
m 


1 - ^ 

a 
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Cr^(t) = - A + (A + C^{0))e 


at 


= i “ ^ + 


A - I + 8e '>3 


“a^t^-t^) 


at. 


= - A + (A + Cj^(0))e 


C (0) = 
m 


- A + 


^ t 
a 


d "a(t4-t3) 

A - I + Se ‘3 


-a{t3-t^) 


-at. 


Each of these equations, with the initial conditions 
substituted into their proper place, yield the complete 
set of mp boundaries shown in figure 2.1. The Horowitz- 
Sidi technique could now be applied to these boundaries 
to complete a design. 

Because the problem is a nonminimum-phase one, 

designing to this set of bounds is not sufficient to 

insure a satisfactory response. A set of frequency based 

magnitude specifications alone is not enough to control 

the response. From the equation describing the response 

due to one RHP zero, C(t) = C (t) - i 6 (t), it is clear 

xii d in 

there must also be bounds on the derivative of the 
minimum— phase response to insure satisfactory nmp res- 
ponses. If both C^(t) and C^(t) fall within their bounds, 
C(t) will also satisfy its specifications. This second 
set of boundaries is easily derived from the boundaries 
already found. It is found by solving this equation 
where all variables are known: 
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Cr^(t) = a[Cj^(t) - C(t)] 

Using the Horowitz and Sidi technique on these two 
sets of tolerances will yield a design for the original 
set of bounds. This is a fairly direct method for de- 
signing to a set of nonminimum-phase specifications as 
long as there is no ignorance in the position of the 
RHP zero. If there is a range of positions where the 
RHP zero could be, designs must be made over the entire 
range and the design that meets all conditions must be 
chosen. This would mean doing many designs before the 
needed one is found. The only time this would be un- 
necessary is when a worst-case set of specifications 
could be designed for and this design would be satis- 
factory for all other possible positions of the RHP 
zero. Since a worst-case would not be obvious, many de- 
signs would probably be required for each set of bounds. 
This is not a practical method of design when there is 
ignorance in the RHP zero. 

2 . 4 Method * s Worth 

This method becomes very impractical when there are 
multiple RHP zeros that must be designed for. More and 
more sets of boundaries are required to specify the al- 
lowed variation of each additional derivative that is 
required as the number of RHP zeros increase. Not only 
do the number of boundaries increase, but the equations 
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that must be solved to find the boundaries become in- 
creasingly complex. Thus, for multiple zeros, the method 
becomes much too complicated and involved to be useful . 
This method holds little promise for a more general solu- 
tion because of the increasing difficulty with more than 


one RHP zero. 



CHAPTER III 


FREQUENCY DOMAIN METHOD 

3.1 Introduction 

This method is to enable a priori design to meet 
the nonminimxim-phase (nmp) tolerances specified. It was 
to be similar to the procedure used to design for 
minimum-phase systems (mp) except both magnitude and 
phase tolerances of the frequency response must be 
studied, but a simplifying discovery was made during 
the investigation. This will be clarified later in this 
section. In a rap system, the magnitude uniquely defines 
the transfer function, but in the nmp system this is 
only true if the position of the RHP zero is known. If 
a design procedure utilizes both the magnitude and 
phase tolerances of the frequency response, several 
designs will not be required as was true with the time 
domain method. 

3.2 Method 

The first step of this method is to transfer the 
nmp step response bounds into magnitude and phase bounds 
in the frequency domain. This is accomplished by approxi- 
mating whatever the final transfer function may be by 

the following; _ 

Oil b 

- (s-a) 

T(s) = 2 -S (3.1) 

(s+b) (s + 2zo) s + w^) 
n n 
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If "a” takes on the values of the RHP zero of the plant 
and z and are allowed to vary over many values, many 
representative step responses can be examined to see 
which ones satisfy the required time domain boundaries. 

A computer program called TEST2 has been written to per- 
form this task. It varies all the parameters of the 
transfer function and finds all the responses which fall 
within the input time domain bounds , The bounds are de- 
scribed by a series of "if statements." All parameter 
sets yielding satisfactory responses and the parameters 
whose responses "hit" the boundaries are printed out. 
These "hits" are stored and when many are found, they 
are used to approximate the time domain bounds by equi- 
valent frequency domain tolerances, magnitude and phase. 
The maxima and minima composite curves and their dif- 
ference are printed. This allowed variation of |AT(ju)| 
is what is required to use the Horowitz and Sidi method 
of design. 

3 . 3 Examination of Method 

To examine the method and to gain insight into 
what is really necessary to complete a design, various 
phases of the general method will be examined. For this 
a useful set of step response boundaries must be found. 
These bounds would be given as specifications in a real 
problem. To accomplish the above, the far pole in 
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equation 3.1, (b) is set at -20, and and z were made 

one for a critically damped system. The RHP zero (a) was 
set at 10 so as to have little effect on the step res- 
ponse. A set of boundaries were than centered about this 
response. The boundary size was then adjusted to include 
some responses and to exclude others. The boundaries had 
to be narrow enough so the boundaries are not trivial 
and cannot be exceeded, but wide enough so they are not 
impossible to meet. The TEST2 program was run for the 
set of bounds shown in figure 3.1 with a RHP zero at 
.75, 1.0, 1.5, 3.0, or 8.0. Figures 3,2, 3.3, and 3.4 
show the output plots of the allowed magnitude variation, 
the allowed phase variation, and the acceptable regions 
of LHP pole variation that yield satisfactory step res- 
ponses . 

The curves for allowed phase variation are explained 
by their region of acceptable pole variations. Since the 
phase variation is entirely due to the pole variation, a 
plot was made for the phase due to the pole locations 
possible in figure 3.4 and is shown in figure 3.5. The 
frequencies where the poles may be on the Bode plot 
account for the behavior of the phase curves. These phase 
curves will not be required in the final design method 
since only one value of the RHP zero will be designed 
for. Since its position is known, magnitude alone is 
sufficient to determine the transfer function. Since the 
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curves do agree with what is expected though, they do 
serve a purpose in that they verify the program’s 
validity. 

The regions of acceptable pole variation (figure 
3.4) get larger as the RHP zero moves further out, only 
up to a certain point, then they start getting smaller. 
The behavior of these curves must be explained. First, 
the lower portions of the bounds will be examined. A 
transfer function with a RHP zero far from the origin 
will have little undershoot and will approach a minimum- 
phase system. One with a RHP zero very close to the 
origin will have a very large undershoot. The allowable 
regions for the complex poles relate to the speed of 
the system responses for various zeros. Clearly, a sys- 
tem with a RHP zero at one will require a faster res- 
ponse than one with the zero at ten in order not to 
violate the lower response boundary (figure 3.6). This 
means a slower system, smaller o)^, will be acceptable 
with the RHP zero further out. This means that as the fhp 
zero moves further from the origin, pole variation is 
allowed at lower frequencies. This explains the lower 
portion of the region. 

The equation relating the mp and nmp responses 
must be examined to explain the behavior of the upper 
boundary. As long as the lower response bound is the 
one that will be violated, C^^^Ct) must be kept small for 
a zero close to the origin in the equation 
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C(t) - C^<t) - i c„(t) . 

If the zero is qIoso to the origin, — is large, there- 
fore to make sure the undershoot does not violate the 
lower bound, the response must be slow to keep C^(t) 
small. If the zero were further out, "a" is larger, 
then C^(t) can be larger without violating the lower 
bound. Therefore, while the lower bound is still the one 
being violated, the allowable region increases in size 
as the RHP zero moves further out, a faster response is 
allowed. 



Figure 3.6 Lower Bound Violation by Nonminimum- 
Phase Responses 
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At some point the response will be so fast, it will 
violate the upper bound. This occurs as the zero goes 
far enough out to allow a response that rises faster 
than the upper boundary does. This has been verified by 
examining various responses to see which portions of the 
boundaries are violated. When this starts to occur, the 
allowed pole variation regions begin to close in. The 
regions begin getting smaller when the upper boundary 
becomes the critical one. Figure 3.7 shows that the res- 
ponse with a RHP zero at one will be satisfactory, where- 
as the response with the zero at ten will be unsatisfac- 
tory. As "a" increases now, the system must be slower to 
prevent violation of the upper boundary. 
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3.4 Design for Worst Case 

Given a set of variation curves as shown in 

figure 3.2, can a design for a certain set of specifica- 
tions and one location of the RHP zero be satisfactory 

as the position of the RHP zero changes? Consider a 

* 

plant P(s) which can be written as P(s) = P (s) (s-a) , 
where P (s) is independent of "a", the RHP zero. In 
using the Horowitz and Sidi design technique, it is 
necessary to have a plant template which can be used in 
connection with the Nichols Chart to design a compensa- 
tion G. To gain insight into the effect of the RHP zero, 
use a plant template for each value of "a". For increas- 
ing values of "a” 


P (s) = P (s) (s-a) 


P(jw) = P (ju) 


^o^+a^ ^180® - tan ^ ^ 


this equation clearly shows that both the magnitude and 
phase of the term (s-a) increases. This means that if a 
plant template is made for each value of "a" the tem- 
plates will be identical but will have greater magnitude 
and phase positions with increasing "a" (figure 3.8). 

The worst case that can be designed for is when . the 
specifications are the seune over the range of RHP varia- 
tion, In this case the same variation must be met whereas 
other times, as the zero increases, more variation is 
allowed (figure 3.2) and the control can be decreased for 
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this case. Certainly the most difficult problem is when 
the strictest set of specifications must be met over 
the entire range of parameter variation. 

If designs for the compensation G(s) were made 
using the plant templates of figure 3.8, the same L(j(o) 
boundaries would result for both zeros if the same set 
of worst-case |AT(jto) | specifications were used. Figure 
1.1 explains the meaning of G(s) and L(s). This happens 
because the size and shape of the plant templates are 
independent of the location of the RHP zero, only the • 
plant template values vary (|P^|,£P^). The Nichols Chart 



Pj^ (s) = P* (s) (s-a^^) 

P 2 (s) ■ P* (s) (a-aj) 
where a^^ < a 2 

Figure 3.8 Template Variation 
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boundary on L{jto) is found by sliding the plant template 


along a constant phase line until it fits between the 
c^^^ves that correspond to the allowed varia- 
tion in T(jw). Since the same template is used for var- 
ious "a's", and the same variation in 

lowed, the same L(j(i>) boundaries will result. The com- 


ious "a's", and the same variation in 


is al- 


pensation G(jw) boundaries are found by: 


iG(jo)) 

'db * 1 

L( joj) 

Idb 

- |P(jto) 1 

db 

|G^(jw) 

Idb “ 1 

L(jw) 

1 db 

- iPj^(jw) 

Idb 

|G2(jw) 

Idb “ 1 

L ( jo)) 

Idb 

- |P2(jw) 

Idb 

from before | 

Pl(jw) 1 

db ^ 

1^2 

l3“) Idb 


which implies 

|Pj^ (j«) 

Ub * 

A = 

- |P2<jia) 1 

where A > 0 db. 

Idb 

+ |Pl(j 

1 db 


iGjfjto) lab 

+ |P2'3“)ldb 

|G 


= 

|G2(jM)lab 


which implies 


'>idb ^ 

> l< 

Idb 



The same argument also will yield /Gj^(jw) > /G^ijdi ) On a 
plot of |G(jci))| vs £G(jaj), this means each corresponding 
point on Gj^{ju)) will be higher and to the right of G 2 ( jw) 
(figure 3.9). If Gj^(ju) is above G 2 (jw) at every point, 
then a design for the smallest value of the RHP zero will 
be satisfactory for all positions of the zero. If G^(jw) 
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is lower than G 2 (ju) at some phase angle, even though 
each corresponding point is higher (figure 3.10), then 
the design must be more complicated because more than 
one location of the RHP zero must be designed for. 

When the design is completed, one set of magnitude 
and phase of L(jai) on the L(jw) boundary is realized. 

The question whether the G(jto) boundaries cross as in 
figure 3.10, is equivalent to asking another question. 
Utilizing the compensation G(juj) that yielded the point 
on the L{jo)) boundary with one position of the RHP zero, 
will the L(jw) found as the zero moves also be on or 
above the L(jti)) boundary? If the Gj^(jo)) used to get the 
point on the boundary is below the G 2 (jw) boundary, 
then the point L 2 found with a 2 and Gj^(joo) is not neces- 
sarily above the L(jw) boundary. 



Figure 3.9 Desirable G Boundaries 
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It is known that if Lj^ is the L(ju) designed for 
with Gj^{jw) and a^^, that the resulting point L 2 as the 
zero moves to a 2 » where 32 > a^^, using the designed 
compensation will be above and to the right of 
on a |L(jto) I vs /L(jw) plot. Whether L 2 is above or be- 
low the L boundary, h 2 ^ or L 22 iii figure 3.11, depends 
on the lead added by (s-a) at that frequency by the in- 
crease in "a" and the slope of the L boundary. For these 
reasons it cannot be proven when it will or will not 
happen, but an argument can be given that it will not 
happen in most problems. 

In the majority of problems specifications will not 

a, 

be met at frequencies higher than w = j , where "a^" is 
the minimum RHP zero location. At frequencies much less 
than almost pure gain is added as the zero moves 

out, eq, 3.2, ^ = 0, This means that at frequencies 

P(s) = P*(s)(s-a) P^(jo))(l - (3.2) 

where w << a, the design for a^^ will certainly be sa- 
tisfactory for 32 * As u) approaches a^, more and more 
lead is added as the RHP zero increases to a 2 . For most 
plants, and an exception has yet to be found, the L 
boundaries appear as in 3.12. 

It is clear that if the phase added by the increase 
in the RHP zero is to yield an L(jo)) that falls below 
the bounds, it would happen at low frequencies where the 
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boundaries are most peaiked. It is at these frequencies 
though where the change is almost pure gain and there is 
no problem. As the frequencies increase^ it becomes 
harder to be above some and fall below the boundary 
also since they flatten out and actually start going ne- 
gative. Many examples have been tried and in all cases 
a design for the smallest zero will be satisfactory as 
the zero varies. This argument at least makes this 



Figure 3.12 Typical L(jw) Boundaries 
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theory more plausible and the examples fail to expose 
possible problem. A rigorous proof is not possible 
since it depends on so many circumstances. Certainly if 
the range of RHP zero variation is small or if the fre- 
quency range requiring feedback is much smaller than the 
smallest RHP zero, a design for the smallest RHP zero 
will be satisfactory for all allowed positions of the 
RHP zero. In the majority of problems this will be true. 
A quick test would be to find all the G(jcjo) compensation 
curves for possible RHP zero locations to see if the 
G(jw) boundaries ever cross. No plant has been found yet 
where this has happened. 

Whatever the design, it must of course meet the 
tightest set of specifications. Thus a design for the 
smallest RHP zero and the lowest set of allowed lAT(ju)| 
specifications will meet the time response specification 
as the RHP zero takes on other values. As shown in 
figure 3,2, the smallest RHP zero also has the lowest 
set of specifications. This is always true from low fre- 
quencies up to some higher frequency where the specifi- 
cations may cross, i.e., a design for as [3, 8], a=3 is 
lowest until w « 3.5. 

3.5 What are the Lowest Specifications 

From figure 3.2, it is clear that the smallest RHP 
zero has the lowest specifications at least at fre- 
quencies where to < j . From studying figure 3.13, a 
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straight-line approximation of the |AT(jtjo)| curves for 
the regions of acceptable pole variation, figure 3.4, 
some general rules may be found. First, the RHP zero, aj^, 
yields the lowest bounds at least until the frequency of 
the maximum position of a varying pole if it is less 
than the value of the RHP zero. If this frequency is 
greater than the RHP zero, the smallest zero, a^^, will 
yield the smallest |AT(jo)) | at least till w = a^^. The 
smallest value of "a" will yield the lowest specifica- 
tions at all frequencies if the allowed region of para- 
meter variation of a 2 encloses that of a^^, where aj^ < a 2 . 
This is shown by ae [.75,1.0] in figures 3.2 and 3.4. 

The specs for a = .75 are lower than those of a = 1 at 
all frequencies, and the poles of a = 1 could vary at 
lower and higher frequencies than those of a =* .75. 

These guidelines may help in deciding which RHP 
zero yields the lowest specs over the frequency range 
of interest. The surest way is to run program TEST2 on 
several positions of the RHP zero and see what specifi- 
cations are lowest. These specifications are required 
in the design anyway. 

For any set of specifications a design for the 
smallest RHP zero will be satisfactory for all loca- 
tions of the RHP zero. A design for the smallest set of 
allowed variation specs on |AT(jw)j will certainly be 
satisfactory when more variation is allowed. Therefore, 
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a design for the smallest RHP zero and the lowest set of 
specs will yield satisfactory step responses as the zero 
is allowed to vary. 

3 . 6 When Cannot the Design Be Completed 

One major drawback of a system with a RHP zero is 
that the bandwidth is limited. This means that an ar- 
bitrary set of sensitivity or disturbance specifications 
are not realizable. Some plant and specification combi- 
nations may require a large bandwidth, but this in turn 
results in an unstable system. Therefore, it would be 
nice to have some guidelines as to when the design can 
or cannot be completed so a lot of time is not wasted 
in trying to find a stable design when it is not possible 
for there to be one. 

One of the steps in the design is to ascertain the 
frequency range over which feedback is required to meet 
the specifications. This is done by seeing where the 
curves of |AT(jw)| = |AP(j(o)|, Up until w^^|AP(jw)| > 
|AT(jo))|, This means feedback is required up to to 
constrain the variation of the transfer function due to 
the plant variation to the specifications required in 
|AT(jo))j, The curve |AT{jw)| is obtained as before by 
the program TEST2 which takes a composite Bode plot of 
the transfer functions which yield satisfactory step 
responses. The curve |ap| is the maximum difference be- 
tween all the Bode plots of the plant as all parameters 
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but the RHP zero vary. 

The frequency where these curves cross roughly cor- 
responds to the bandwidth of the system. If this fre- 
quency is less than half the value of the smallest RHP 

zero, the design can be completed satisfactorily. If 

a, 

this frequency range is ^ < w < 2a^, where "aj^” is the 
smallest RHP zero, the design may be completed, but it 
will be very difficult to do if it can indeed be done. 

If the crossover is at a frequency greater than twice 
the smallest RHP zero, the design is impossible. This 
should give a rough idea of how difficult a design is, 
if it is possible. 

Another means of judging if the design can be com- 
pleted concerns the Nichols Chart, The loop transmission 
can be written as: 

L(s) = L^(s) (s-a) = L^(s) (s-a) jlJll = L^(s) {s+a) [|^] 
L(s) = L^(s)A(s) 

where A(s) is an all-pass function. 

figure 3.12 was the Nichols Chart boundaries on L(s) , 
a plot of those for would be shifted to the right 

by the phase contributed by the all-pass function at each 
frequency. If there is a disturbance in the system, it 
must be damped out. The damping can be approximated by 
a dominant pair of complex poles whose damping factor z 
can be related to the disturbance overshoot allowed. 
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This can then be related to a magnitude peaking curve 
which yields a forbidden region on the Nichols Chart. 

The region specifying the damping of a disturbance, the 
ellipse-like shape in figure 3.12, will move also, but 
it will be distorted because each point on it corres- 
ponds to a different frequency. This is clear if the L 
found wraps right around this region. Each point will 
move a different amount because the all-pass contributes 
a different amount of phase at each frequency. 

If the last boundary corresponds to a frequency 
cj > 2a^^, the damping region will probably cross the zero 
degree line, figure 3.14. When this happens, a design 



I 
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cannot be made because it is not possible to come down 

in magnitude any distance in decibels with phase lead. 

^1 

If 0 )^ < ^ , it is clear a design can be found. If 
®1 

2a^^ > , it will be more difficult to stay out- 

side the disturbance region. This method of looking at 
the problem yields the same results as before but by 
being concerned with the damping of a disturbance signal. 

If a lot of gain variation is present in the plant, 
the last specification curve must be at << a^, because 
it will take a wide frequency range to clear the dis- 
turbance region. This means that if was near a^^, the 
disturbance region would almost certainly cross the zero 
degree line, and probably over a large magnitude range. 
This would definitely make the design impossible to com- 
plete (figure 3.15). 
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3 . 7 Siinunary of Design Procedure 

This section will present the design steps of the 
proposed procedure. First, find the set of frequency 
domain specifications on |AT(jw) j by running a program 
like TEST2 which approximates T(s) by a third order 
function and searches for valid transfer functions. 

Next, ascertain the set of specifications that will be 
designed to, the zero designed for, and the frequency 
where feedback is no longer required. This is usually 
the smallest set of tolerances, the smallest possible 
RHP zero, and the frequency where I AP(jWj_) [ = lAT (jtOj^) I . 
With this information use the Horowitz and Sidi design 
technique to find the boundaries on the loop transmis- 
sion up to 

Now by some technique, possibly trial and error, 
find an L(jw) that meets these requirements and yields 
the desired amount of disturbance damping to the ac- 
curacy desired. Knowing the plant and the loop transmis- 
sion, the compensation G(s) can be solved for. With this 
G(s) and the varying plant a set of actual |AT{jto)| 
specs can be found. The prefilter F{s) can be solved for 
which shapes the transfer functions to the same shape 
and bandwidth of the specifications. This completes the 
design and should yield satisfactory step responses. 
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CHAPTER IV 
FIRST EXAMPLE 


Now that a design technique has been postulated, 
it is necessary to work a meaningful example that will 
test this postulate. The first step in postulating an 
excunple is to have a valid set of time domain specifi- 
cations. The specifications must be tight enough to make 
the example meaningful, non-trivial, and yet not too 
tight so as to make it impossible to complete the design. 
This was accomplished by approximating the transfer 



as was done in the earlier theory section. This approxi- 
mates the transfer function essentially by the RHP zero 
and a dominant pair of complex poles. As before, the 
time domain bounds are centered around a critically 
damped response of T(s) and then situating the bounds 
so some responses are confined to the interior and other 
responses exceed the bounds. The resulting bounds are 
shown in figure 4.1. 

To use the Horowitz and Sidi design technique a 
set of I i^T I specifications in the frequency domain 
are required. A program, TEST2, was written to vary z 
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and in T(s) and find many T(s)'s whose step responses 
satisfy the time domain specification. Frequency domain 
plots of those that "hit" the boundaries yield the equi- 
valent frequency domain bounds of the original time 
domain specifications. The required |ATj^(jto)| specifica- 
tions are the differences between the minimiam and maxi- 
mum plots at each frequency. 

Specifications were obtained for several RHP zeros, 
figure 3.2. Since the allowed |AT(jw) j variation is al- 
most the same for ae[1.5,8], it was decided to let the 
RHP zero vary from one to eight so the theory would be 
tested to the fullest by designing over a wide range of 
allowed variation. That is, if the theory were wrong, it 
should certainly fail when there is a large difference 
in specifications over the allowed range of RHP varia- 
tion. 

There is a rule of thumb that the maximum obtain- 
able bandwidth in a nonminimum-phase system is (j = j , 
where "a" is the position of the RHP zero. If the RHP 
zero varied from .75 to 8. instead of 1 to 8, this would 
mean feedback would be required over the range 
we[0,.375] and the tolerances for a = .75 allowed very 
little plant variation. The remaining parameters of the 
plant could not have any meaningful variation without 
requiring a larger bandwidth. For these reasons a = 1 
was chosen as the minimum RHP zero for this problem. 
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By examining the allowed j AT j variation a plant 
was chosen which would not require feedback over tij = .25 
or ^ . The plant chosen was: 


= ififbT ' 


ae [1, 8] 
be[.015,.15] 


Feedback is not required once I AP I , is below I AT I 

' 'a=l ' ‘a=l 

where figure 4.2. Since the 

variation in the overall transfer function is due to the 


variation in the plant, once the plant variation is be- 
low the allowed variation of the transfer function no 
further control over the plant is required to hold it 
there. Now the example has been completely defined con- 
sisting of a set of time domain boundaries and a plant 
with RHP zero and LHP pole variations. The remainder of 
this example will show how the design technique is ap- 
plied to this type of problem and how well it solves 
the problem. 



Figure 4.2 Range of m Requiring Feedback 


log CO 
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The first step in using the Horowitz-Sidi design 
technique is to arrive at a set of boundaries on the 
loop transmission Lj^(j(o) where Pj^(jo)) is the plant at a 
set of parameters yielding maximum phase lag at all fre- 
quencies. If the compensation G that is used with this 
yields an which does not violate the disturbance 
damping region and is stable, then the plots of the 
other L's resulting from all other sets of plants para- 
meters will be to the right of on the Nichols Chart. 
This means that for all plant conditions, the resulting 
L is stable and has at least the required damping of a 
disturbance. This is why this P^^ is used in the design. 
In this problem /P(jw) can be expressed as: 

/P(jw) = 90° - tan"^ ^ - tan"^ ^ . 

It is clear that the most phase lag is available when 
both a and b are at their minimum conditions, i.e., 
a = 1/ h = .015. The design technique indicates that a 
design for the smallest value of the RHP zero will work 
as the zero varies. Therefore the smallest value of "a" 
will be the only value of concern throughout the design. 

A program implementing the Horowitz and Sidi tech- 
nique was used to yield a set of boundaries that 
meet the input |AT(jw)| specifications by using a plant 
template that describes the allowed plant variation. 
(Figure 4.3) Any compensation which in connection with 
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yield a set of magnitude and phase values of 
L(jw) that fall above these boundaries at all frequen- 
cies will more than meet the allowed |AT(juj) | varia- 
tions, The L(jw) meeting all the boundaries is found by 
trying various pole zero combinations with the RHP zero 
at its lowest value until an acceptable solution is 
found. 

Also shown in figure 4.3 is a small elliptical 
shaped region near the -180“, Odb point. If there is a 
disturbance in the system, it must be damped out. If 
the damping can be approximated by that of a pair of 
dominant complex poles, the damping factor z can be re- 
lated to the disturbance overshoot that can be allowed. 
This z can then be related to a magnitude peaking curve 
which shows up as a region on the Nichols Chart. The 
region used here is approximately a 3db curve and cor- 
responds to an allowable 25% of disturbance overshoot. 

The optimum L(jco) would fall on each L-boundary at 
each frequency and follow the 3db boundary without 
entering it. Using the trial and error method of design, 
this is the most time consuming portion of the design. 

As the frequency approaches that of the RHP zero, phase 
lag is introduced with gain. This makes it extremely 
difficult to stay below the 3db region. If feedback is 
required at half or three-fourths the value of the RHP 
zero, it may be impossible to find an L(jw) that acts 
in the desired manner. 
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The L^(s) for this example is 


-640000(5-1) (s+,15) (S+.Q15) (s+.4) 
s(s+.015) (s+100) ^ (S+.02) (S+.24) 


Li2(s) 


fe^+.22s4-,0625) (s^+1.4s+l) (s^+. 105s+ ♦ 005625) 
(s^+.35s+.065) (s^+1.85+1) (s^+. 135s+. 005625) 


Ij^(s) — ( s ) ( s) 

Once a suitable Lj^(s) is found, it is necessary to 
solve for the required compensation G(s). This is done 
as follows: 

L. (S) 

L^(s) = G(s)P^(s) => G(s) = p 

For this problem, the following G(s) resulted: 

-640000. (s+,15) (S+.015) (s+.4) 

(s+100)^(s+.02) (S+.24) 

(s^-i-.22s+.0625) (s^-t-1.4s+l) (s^-h. 105s+. 005625) 
(s^+.35s+.065) (s^+1.8s+l) (s^+,135s+. 005625) 




G22(s) - 


G(s) = (s)G 22 (s) 

If this were the optimum G(s) , yielding an L(jco) 
on the boundaries at all frequencies, it would yield a 
|AT(jo))j, as the parameters of P varied, which would be 
identical to the |AT(ju>)| specified. Since all points of 
L(jw) are on or above the boundaries, the resulting 
|AT(j(o)| for the example, call it |AT 2 f, should be the 
same or less than that specified, called | AT^ | . Figure 
4.4 compares these two curves. It is clear that the 
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curves are not as expected. The error is not in the 
method of designing for the smallest zero since the zero 
has not been allowed to take on any other values . After 
examination of the program which finds the roots of T{s) 
from L(s), an error is clear. This program takes the 
^oots of L(s) and outputs the poles of T(s) . For one run 
there is a pole and a zero of L at .015 and therefore 
there should be a pole and a zero of T(s) at the same 
point. The program says the pole is at .0155 instead of 
.015. This error alone leads to a |AT 2 | that is in error 
by easily the amount encountered. The method used by the 
subroutines of this program lead to a significant error 
when poles are used at such low frequencies. Since there 
is an error, but an unavoidable one for this example, 
the example will continue. This error does not appear 
to be a significant one. 

ftnal step in the desxgn is to find the prefil** 
ter that will put the variation in the shape and band- 
width of the allowed region indicated by the components 
of I I . The compensation G holds the plant variation 
to an acceptable value and F shapes it to match the T 
specifications* Firom the equation below ^ a way of cal- 
culatinq F(s) is to use the maximum curve making up 
IaTj^I and subtracting the maximum curve of | AT_ | in 
decibels. This difference can be realized by some P(s) 
that will make the curves coincide. If the optimum G(s) 
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had been obtained, 




F(jw) I = 


Tj^(ju)) 


|F(jc«) = |T^(jo.) - |T2<j») 


db 


this F (s) would force the minimum curves of and 
to be the same also since |ATj^| = IAT 2 I in this ideal 
case. Since the G(s) used is only close to ideal, two 
prefilters would result, but they are so close to being 
the same that either should yield a satisfactory res- 
ponse. Figure 4.5 shows these curves and the F(s) 
realized which was; 

F(s) = 53.352 (s^ + .26s + .0625) 

(s+2) (s+3) ^ (s^ + .35s + .0625) 

Now the design for the closed loop structure of figure 
1.1 is complete. Figure 4.6 shows the step responses of 
the system under various plant conditions. It is clear 
the design worked very well. The allowed region varia- 
tion is utilized to the fullest and all responses are 
acceptable. 

Some things should be pointed out about the example. 
Usually the simplest problem is gain variation, but this 
was not possible here. The specifications on |AT^1 were 
so low in this frequency range that only an extremely 








53 


low gain variation could be used. Since it would be so 
impractically small ^ it was left out of this example. 
Another example with larger tolerances must be worked 
in order to demonstrate the technique's validity in 
problems with gain variation. 

This example was to design over a range of RHP 
zeros where the |AT{jw) | specifications vary greatly. 
Variation from one to eight demonstrates how correct the 
premise is of designing for the smallest zero. If the 
specifications were RHP zeros from three to eight where 
the tolerances are almost constant, there could still be 
doubt in the usefulness of the theory here extended. A 
range of specifications will not be of main concern in 
the next example since this one showed how effective 
the design technique is. 

Another example will therefore be worked in this 
area to investigate gain variation. A higher frequency 
range will be worked with to diminish the effect of 
errors encountered in the subroutines of the program 
that finds the roots of T(s). 

This should help clarify whether the subroutines 
are the cause of the error or not. Also another example 
will stress the effectiveness of the method. Modeling 
time domain boundaries by a third order transfer func- 
tion gave a very adequate set of frequency specifica- 
tions on I AT ( jo)) I . 



CHAPTER V 


SECOND EXAMPLE 


Whereas the first example was concerned purely with 
verifying the design theory be designing over a range of 
specification variation, this example was to answer other 
questions. The first example, although verifying the de- 
sign technique, had raised some questions that needed 
more than a possible explanation. First, was the exces- 
sive variation of the designed transfer function ( | AT 2 | ) 
really caused by the inaccuracy of the subroutines of 
the root finding program? Second, gain variation was im- 
possible due to the low specs that had to be met in the 
first example. Does the design method work for gain 
variation, or are further problems encountered? To an- 
swer these questions, this example has gain variation 
and deals with a higher frequency range to help decrease 
the error in the subroutines. 

This example uses the same time domain specifica- 
tions as the first example, but the plant is changed to; 

ae [3,8] 

" sjs+b) ' Ke[l,2.2] 

be[.075, .75] 

This plant needs feedback to meet the design specs up to 
a frequency of .72 RPS. This was found by seeing where 
the specs become lower than allowed 
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specs (Figure 5,1). 

The next step is to find the L-boundaries for the 
plant over this frequency range. These are found by the 
program implementing the Horowitz and Sidi method and 
are shown in figure 5.2. Also shown in the figure is 
the region that must be avoided in order to permit only 
25% overshoot of a disturbance signal. 

Finding a loop transmission L(s) that meets these 
boundaries is the most difficult part of the design. 

The optimal solution is to meet each boundary at each 
frequency and to come as close as possible to the for- 
bidden region. Shown in figure 4.2 is a plot of the 





2.2(j(o+.75) 

l(jaj+.75) 


Figure 5.1 Frequency Range Requiring Feedback 
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Using this compensation and allowing the plant parameters 
to vary throughout their entire range, holding a = 3, 
will yield the actual variaton of the transfer function 
lATzIdb* example one, the maximum component of 

I AT 2 I will be required to complete the design by solving 
for the required prefilter F(s). Figure 5.3 compares 
I AT 2 1 and the specs of | AT^^ | . If the found in figure 
5.2 had points on all the L-boundaries instead of some 
above and some below, the plots of | AT^^ | and 
I AT 2 I would be identical. When a point is above the 
L-boundary, the IAT 2 I at that frequency will be below 
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that for IaTj^I and the opposite is also true. This was 
the problem encountered in example one and thought to be 
caused by errors due to the techniques used in the sub- 
routines of root finding program. Note that the problem 
is not encountered here as was hoped. Poles do not occur 
at the low frequencies encountered in the first example, 
and therefore if the subroutines were the cause of the 
error they should not and do not lead to erroneous re- 
sults here. This answers or at least supports the sus- 


pected answer to the question concerning the origin of 


the error encountered in the first example. 

Figure 5.4 shows four things, the maximum component 

of l^*^lldb ^"^Im^ ' maximum component of lAT 2 |^j^ ^^2m^ ' 
the prefilter required and the prefilter realized 
F^. As before, Fj^ is simply the point by point subtrac- 
tion of T 2 jj^ from T^^^ as explained by: 
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The prefilter found is; 

F(s) = 3.49(s+2.5)^(s^+.42s+.176) 
(s+.7)^(s+4)^(s^+.98s+.49) 

Thus the design is now complete. Figure 5.5 shows exact- 
ly what is meant by a prefilter F and compensation G 
as to their part in the closed loop structure. To test 
the design, all parameters were varied over their al- 
lowed range and plots of their step responses were made. 
Figure 5.6 shows that again the design technique works 
perfectly. The allowed time response region is utilized 
to its fullest as the economical design should. 

This example accomplished two things. First, it 
again verified how well the design technique functions 
when applied to a problem. Second, it helped answer the 
questions that prompted a second example. The inaccuracy 



T(s) 


C(s) 

RTiT 


L(s) = G<s)P(s) 


Figure 5.5 Closed Loop Structure 













of the subroutines of the root finding program at very 
low frequencies certainly seems to be the cause of ex- 
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cessively large variation of | AT2 [ at certain frequen- 
cies in the first example. The design technique works 
for gain variation without introducing any complications. 

The example again strenghtened the feeling that a 

a . 

design for a system requiring feedback up to w = 

is not an easy problem, but is a workable one. A problem 

a . 

requiring feedback up to — ^ will be much more difficult 
if the trial and error method is used to find the L 
meeting the required L-boundaries . This is the only dif- 
ficult part of the design and requires time in exponen- 
tial proportion to the accuracy desired. 



CHAPTER VI 
CONCLUSIONS 

The first design method presented in this paper is 
not very useful in the majority of real world problems 
because of its inability to handle variation of the 
right-half-plane zero. Its technique may yield a more 
general method or at least give some insight in the 
problems encountered in a nonminimum-phase problem. 

The second method is very effective in presenting 
an a priori design technique for a plant with one right- 
half-plane zero and the plant is not a function of the 
right-half -plane zero in any other way. Necessary exten- 
sions of this method are to handle multiple right-half- 
plane zeros and plants that depend on the right-half- 
plane zeros in more complicated ways than just having a 
zero there, i.e., gain is a function of the position of 
the zeros. 

Some of the simplicity of this method may be lost 
when these extensions are made, but this method allows 
an insight into phases of the design which point out 
short cuts. For example, it made it clear that under 
certain circumstances a design for one position of the 
zero will work for all variations. It also facilitated 
the presentation of criteria explaining when a design is 
not possible by this method. It seems very likely that a 
pursuit of the nonmmimum— phase problem along these 
lines may yield a general design technique. 
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Appendix A 

DIGITAL COMPUTER IMPLEMENTATION OF TIME DOMAIN TO 
FREQUENCY DOMAIN TRANSFORMATION 

A. 1 General Description 

The program which yields the |At| specs necessary 
to make the design and the maxima and minima curves of 
T to calculate a prefilter F is presented here. It is 
coded in FORTRAN IV source language and has been exe- 
cuted on a Control Data Corporation 6400 computer. A 
series of IF STATEMENTS in the subroutine INVLAP de- 
scribe the step response boundaries. A "hit" is a res- 
ponse that at some time came to within ±.01 of the 
boundary without ever exceeding this tolerance. This 
tol®i'3nce and the amount the complex pair of poles are 
incremented each time can be adjusted to suit the ac- 
curacy desired. 

The program varies the complex poles of 

( 0 ^ 20 
(s-a) (-ii- — ) 

T(s) = —2 ^ 

(s + 2zco^s + 0 )^) (s+20) 

^tnds the region of pole variation that yields step 
responses within the boundaries. It will print out the 
first and last set of parameters giving satisfactory 
responses for a given real value of the poles, varying 
the imaginary part. The sets of parameters that yield 
"hits" are also printed out. Finally the 


program prints 
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out the maximuin, mininmin, and difference of the magni~ 
tude and phase bode plots of all the hits. All this in- 
formation provides the insight necessary to design for 
one set of parameters and yield satisfactory responses 
as the parameters vary. 

The following is a logic flow chart of the program 
and a complete listing of the program with a sample set 
of output. 
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A. 2 Flow Chart of Time to Frequency Domain Transformation 


(START) 


INITIALIZE CONSTANTS 


READ ZERO - RHP 


INITIALIZE REAL PART OF COMPLEX POLE PAIlH 


INTIALIZE IMAGINARY PART OP COMPLEX 
POLE PAIR 

INCREMENT REAL PART BY .1 


Print Last, 



INCREMENT IMAGINARY PART BY .05 
CALCULATE HIGH FREQUENCY GAIN 
CALCULATE TIME RESPONSE 


IS \ 
STEP RESPONSi 


_X^IS FIRSTS. 
IN RESPONSE FO^ 

\this real 

\ PART / 


FIRST 

NO y/^ME WITH THiJ 
“X. REAL PART 


IMAGINARY 
s^RT > 15 
\ 7 / 



PRINT OVER 15 
FIND MIN & MAX OF 
BODE PLOTS OF ALL HITS 


PRINT MIN, MAX, AND A 
OF MAGNITUDE & PHASE 
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A, 3 



PR09RAM TEST2(InPUTiOUTPUT) 

dimension PD( 20 )fRNI 20 )»A< 10 l« 9 )tCN( 20 t*P( 5 )fZ(S>tPCl 2 ( 10 )« 
j 2 B( 10 It 9 ) •DAT( 200 «A) >ZC( 2 « 10 ) 

COMPLEX RNiRO 
! 10 ITEMaZ 

I CN( 1 )> 0 * 

IPINl -0 

Ad.Da.I 

ROa>ACMPLX(0.0,0*9> i 
RD(2)acMPUX(*ZO.OtO.O) 

20 READ 900 *ZERO 

j 900 FORMAT(F 10 » 6 ) 

’ PRINT 901 (ZERO 

’ 901 FORMATdXiFlO. 6 ) 

j IHIT-l 

IF(ZERo »LE* .11 00 TO SOO 
IFINlalFINI«l 

I ' RN(l)aCMPLX(ZCR 0 i 0 .e) 

i POLERa .3 

30 POLEIaO. 

poler«poler«.i 

FRST- 0 . 

35 P0t.EI*P0LEI*.05 
iNaO 

RD< 3 »bcmplX(-P 0 lER»*P 0 LE 1 » 

RO<*>acMPUX(-POLER»POLEl> 

^ SFa •l 30 ./ZERO>*(POt.ER** 2 *POCEI*aZ) 

lFlRET. 3 .) 35 f 30.40 
♦0 ITEMlaxMIT -1 
Ob 39 iBidOO 
A(|» 2 )a*lOOO* 

A(Ii 3 )al 000 . 

B(IiZ)a- 1000 * 

39 S(l* 3 )al 000 . 

00 100 lalilTEMi 

P< 1 )» 0 AT(I.H 

Zll)aOAT(l> 2 ) 

PC(ltl)aoATlI. 3 ) 

PC 12*1 1 boat (I. 4} 

bal. 

NPal 

NPCal 

NZal 

OELa.l 

NZC *0 

Sao. 

CALL BODE<S»OiNp,P*NPC*PCiNZ. 2 |NZC,ZCoD 9 ia*OELoA,B) 

100 CONTINUE 

1N0W-IFINI«3 

DO 200 lal.lOO 

A(I,lN 0 W)aA|I» 2 }.A|I, 3 ) 

d(l«lN 0 W) 8 B(Ii 2 ).H(it 3 ) 

i 200 CONTINUE 

600 PRINT 910 ' ^ 

DO 300 laltlOO 

PRINT 905 *(A(I*j)tJal* 4 )«<eiI*J)lJa 2 *A) 

300 CONTINUE 

905 F0RMATtlX*7F13*4> 

910 FORMAT(lHl) 

PRINT 910 
90 TO 10 
500 CONTINUE 
990 end 
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Z^&HUf A • ITcHtZHITt *Nf FWST ) 

comment JOUJINC evaluates the INVpSE LAPLACE TRANSFObM Of A RATIONAL 

c function at nt points at intervals ot» 

C INPUT FORMS 

C MS* l« J wiSIliilir:: *! <NN) » > / < ($.RD U > > #. . ,# JS-RO <N0M > 

c '<*2t Sf*CS»*NN*C(2l«S#«(NN-n*t*.^C(NNI*5*C{NN*l) »/ 

® . ns-RO(i) j*.,.*(s«RO(NO)n . 

. INTEOER hit 
HIT»0 

^EfM^FlolA.i^^CTli^aeiE^S^ RESPONSE OUT OF LIMITS AT 

DO T laliNO " * * ’ 

RESIDaSF 

if(nn*eo.o) 00 to 6 

00 T0J2*4)»K 

E 00 3 JbIiNN , 

i RES(I)aRES(I>*(R|)(I).RN(jn 

00 TO « 

A RES(I>a(H0(I)*CN(2n 
00 5 UbEiNN 

5 RES(I|bRES(I>*RocI|*CN(J^1) 

REStnaSr*RES(I| 

6 DO 7 jal.NO 


if(Ueq«j} 00 To r 
RESII)bRCS(1I/IR0II)>R0(J) I ' 

' 7 CONTINUE 

TIbAJI.U 
00 9 taUNT 
CTliO* 

DOS UaliNO 

moiJi ) 1 . 

print i»TlfCTl 

^ A(IiITCH)bREAUcTU 
OS IFCTI «0E« 9*07) 00 TO 790 

^**%iJof))OO^TO*700 *‘***^^”‘ .LE.U.ID.AND. A II, ITEM) , OB. 

102 IFITI .OC* r.) oO TO 125 

^”aOO*TO *®^*‘**’ ••'‘°»*‘*'*TEM) *LE. (1.21). ANO.Ad, ITEM) .OE, 1.64) ) 

1.J ;f,Ti M. ... ,0 TO >» 

112 IF(T1 .OE. 3.) 90 TO 129 

'•••'■•.on .u.iitw 

*^<Tl ,0E. 1.) 00 TO 125 

* ®zo5 To‘j5r^”’ *‘***TEM) ,LE. (.4*T1*.0D) 

125 HIT«0 

IF (FRST *L E. «5 ) 00 TO 130 

127 ONEOAasoRT(POUCR«*2*POLEl**2) ” ' ■ ■ ' ^ 

ZETAaPoLCR/OMEOA 

S?£2V^‘*^°‘-^**»'’®'-E**2ET4,0ME0A 
FAST^O* ' ■ ’ \ ■ 

■ «ET*3* . - . ... ; . , , V . . jo/.; 

00 TO IS 

J2® ll ‘•*0LEI *0E. is.) 00 to 190 - 

Mi ll ii** ** ®® TO 370 

145 60 TO 10 - ■ , * 

ISO RET*4, 

PRINT 073 ■— •- 'X . 


print l.Tl.CTl 



31$ $0 TO f 
IBO INM 
•0 TO 0 
3T0 XN»0 

00 TO 10 

400 IF (AdtITEM) .),£• 1 .31667«T1>.60647) .OR. A<I»ITEM) .OE. 

2|.4«n..0t>> HI7>4 
410 00 TO 310 

SOO If: (AlItlTEN) .lE« (.31667«T1-.60667) *0R. AII*1TEH> .OEi (I.IOH 
2H1TM 

SOS 00 TO 310 - - , 

400 IF (A(idTEM> .lE* (*06) *0R« A(I«1TEM) .OE* a*19)> Htt>4 
OOS 00 TO 310 

too XF(A(I.ITEM) *LE. t.91) .OR. A|I«ITEN) .OE. (1.09)) HIT*4" 

70S 00 TO 310 

7S0 IF (FRST .OE. (.S)) 00 TO 775 
0«e0A"S0R7 (P0LBR**0*F0LBX*»2) 

F ZETA>ROLER/OMEdA 

PRINT 970tPOLER,POLEI*ZETAfOMEOA 

970 format (1XiSHREAu«*F 10. 4»3XfSHlMAaa.FlO. 6. 3XtSHZETAa«F10.0t3Xi 
24HOMESAa.F10.4t4X»5HFIRST) 

FRST»l. 

778 IF (MIT .LE. 2) 00 TO 10 
760 OMEOAaSORT (P0LER*«2*P0LEI*«2) 

ZETAaPOUER/OMEOA 

PRINT 972.P0LERtP0LEItZCTAt0ME0A 

972 F0RMAT|lX»SHHEAU*tF10.6.3XtBHlMAOaiFlO.4»3X«5HZETA>*Fl0.6»3Xf 
26HOMEOAB«F10.4f4X*3HHlTi/) 

MITaO 

0AT(1H1T*1)>20. 

OAT(IMIT»2) — (ZEBO) 

0AT(1M1T,3)«ZET4 

0Ar(IMIT.4)a0MEaA 

IMIT-IHIT*! 

9 TlaTl* 0 T , 

10 RETa2. 

15 HCTUMN 
- END 


SUBROUTINE BOOE(S«0*NPiP*NPCtPC«NZ*Z»NZC*ZC>Oa»Q«OEL)A»e> 

COMMENT SUOROUTINE EVALUATES RATIONAL FUNCTION AT S. 

comment POLEsaP(I). ZeROSaZU). complex POLESaZ£TAt0MEOAapC(l«I)tPC(2fl) 
comment complex ZEROS«ZCTAfOMEOAa2C(liI>.ZC{2tli* ERROR CONSTANTaO. 
DIMENSION P(l)iZ(l)(PC(2»10)iZC(2ilO)»A(10l«3)*S(101i3) 
PIa3*UlS927 

00 100 IBIOal.loO - . 

SbS«DCL - 

AdBIO.DaS 
e(IBI0.1)a$ 

' ■ DBaO**2 

= QaO. 

if(np*lc.o) 00 to 3 

00 2 lalfNP 

, xf(P(I).eb«o.) 00 to 1 

OBaDBBp ( I ) a«2/ ( s#»2«p I I ) «#2) 

Uae.ATANIS/Pll)| 

• 00 TO 2 - - 

rv 1 DSaOB/saa? 

■' OaQ-Pl/2. 

2 CONTINUE 

3 IF(NZ.LE.O) 00 to 6 

j ' " 00 S laliNZ ' . . ■ - • - 

r" IF(Z(I).EO.O.) 00 TO 4 

' 0faD6a(Sa*2»Z(||*a2)/Z(I)«a2 

OaO*ATAN(S/Z(I>) ' - 

00 TO S ' ■ • 

4 OBaOBBSaaZ 

QpO*PI/ 2« ■ ■ t 

5 CONTINUE 

:• .0 iF(npc.le.o) oa to 9 
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I 


J 

i 


DO 8 ialtNPC 

08>08*PC (2*1) •♦♦/ MPC ( 2* I ) ••2-S**2) (PC ( 2* 1 ) •PC ( 1 * I ) *5 ) ••B) 

Q-0-ATAN(2.#PC(l.Il»PC{2,n*S/(PC(2,I)**2-S**2» 

lF(S.LE.PC(2»n> 90 TO 8 * *» c a «i» 

IP«PC(l,I).LT.O.) 90 TO 7 
^ Q«Q«P1 
90 TO 8 

7 9aQ*PI 

8 CONTINUE 

9 IP(N2C.UE*0) 90 TO 12 
90 li I«1*N2C 

0B«08* I (ZC (2*1 ) ••2«S**2)**2*4«* (2C (1*1)*ZC(2«I> *S) P#?) /ZC(2> I>*#4 

a-Q*AT*N( 2 .»ZCIi,n*ZC( 2 ,n*S/( 2 C( 2 ;^#»|-s;i 2 )* *’'*‘^‘*’*^ * 

1^{S.(.E.ZC(2*1)) 90 TO 11 
IP(ZC(1iI>.lT*0.) 90 TO 10 
a»Q*Pi 
GO TO 11 
10 ObQ«P1 


11 CONTINUE 

12 0BBS0RT(AeS(D8l ) 

IP (DB .9T« 0») 90 TO 13 
PRINT 937ilBI0 
937 FORMATdOXillO*/) 

PRINT 901*P(l)*2(l)tPC(l«l),PC(2>l) 

901 FORNAT(1X**P10.6> 

13 Ra0B*C0S(0) 

AMaD8*SIN(OI 

DeB20.*AL0010(08) 

8bO«IB0./P1 

1F(A 11810*2) •9e« 08) 90 TO 50 
4S A(lBlDf2)B0B 

50 1F|A(1b 10*3) .tE» DB) 90 TO 60 

59 A(1BIOi3)bOB 

60 1P(B(1B10*2) *9C* 9) 90 TO 70 
65 B(lBlDt2)BQ 

. 70 IP (B(lB10*3) *|.E* Q) 90 TO 80 
75 B(1B10»3 )b0 
80 CONTINUE 
100 CONTINUE 
RETURN 

_ End 




5*000000 


REALb 

*400000 

IMAOa 

REALb 

•400000 

IHAOa 

REALb 

*400000 

1MA0« 

REAC" 

*900000 

XNAQS 

REAUb 

•900000 

IMAOm 

REAL* 

•500000 

IHA9* 

Real* 

•900000 

IHAd* 

REAL" 

•600000 

JMAG« 

REAU* 

•600000 

1HAG9 

real* 

•600000 

1NAG« 

real* 

•600000 

IMAG» 

REAL* 

■700000 

XHAOtt 

real* 

•700000 

XHA9« 


•350000 

ZETAb 

•752577 

•750000 

ZETab 

•470588 

•600000 

ZETAb 

•447214 

•650000 

ZETab 

•894427 

•650000 

ZETab 

.894427 

•950000 

ZETAB 

.465746 

^•oooooo 

ZETab 

•447214 

•OSOoOO 

ZETAB 

•996546 

•050000 

ZETAB 

•996546 

•lOOoOO 

ZETab 

•986394 

l*000o00 

ZETab 

•514496 

•050QOO 

ZETab 

*997459 

^•OOOOOO 

ZETA* 

•973462 


OMEOAb 

OMEGAa 

.531507 

•850000 

FIRST 

HIT 

OMEGAb 

.894427 

LAST 

OMEOAb 

OMEOAb 

•559017 

•559017 

FIRST 

HIT 

OMEOAb 

1*073546 

HIT 

OMEOAb 

1.118034 

LAST 

OMEOAb 

OMEOAb 

•602080 

•602080 

FIRST 

HIT 

OMEBAb 

•608276 

MIT 

OMEOAb. 

1*166100 

LAST 

OMEOAb 

OMCGAb 

•701783 

1.220656 

FIRST 

LAST 
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REAL* 

REAL* 

tBOOOOO 

tIOOOOO 

xhag« 

XMAO« 

.OSOoOO 

l.OOOOOO 

2ETA» 

ZETA« 

•99B053 

.624695 

OKEOA" 

OREOA" 

.801561 

1.280625 

FIRST 

CAST 

REAL" 

REAL" 

•900000 

•900000 

lMAo« 

IMAO« 

.050000 

.950000 

ZETA* 

ZETA- 

.998460 

.687745 

OMEOA" 

0ME8A" 

.901388 

1.308625 

FIRST 

HIT 

REAL- 

.900000 

IMAQ9 

l.OOOOOO 

zeta* 

.668965 

OMEOA" 

1.345362 

LAST 

REAL" 

REAL" 

l.OOOOOO 

1*000000 

XMA3« 

xnag« 

.OSOoOO 

•950000 

ZETA* 

ZETA« 

.998752 

.724999 

OMEOA" 

OMEOA" 

1.001249 

1.379311 

FIRST 

HU 

REAL" 

l.OOOOOO 

xhag» 

1.000000 

zeta» 

.707107 

OMEOA" 

1.414214 

LAST 

REAL- 

REAL" 

1.100000 

1.100000 

XMAG- 

XHAO* 

•OSOoOO 

.900000 

zeta* 

Z£TA« 

*998969 

.773957 

OMEOA" 

OMEOA" 

1.101136 
i. 421267 

FIRST 

HIT 

REAL" 

1.100000 

XHAO* 

.950000 

ZETA* 

•756823 

OMEOA" 

1.453444 

LAST 

REAL" 

REAL" 

l.EOOOOO 

l.EOOOOO 

|MAG> 

imag* 

.OSOoOO 

.650000 

ZETA« 

zeta* 

•999133 

.816024 

OMEOA" 

OMEOA" 

1.201041 
i. 470544 

FIRST 

HIT 

REAL" 

l.EOOOOO 

XNAG« 

.900o00 

ZETA* 

•800000 

OMEOA" 

1.500000 

LAST 

REAL" 

REAL" 

1.300000 

1.300000 

XHAG« 

imag« 

.OSOoOO 

•OOOoOO 

zeta« 

ZETA* 

.999261 

.851658 

OMEOA" 

OMEOA" 

1.300961 

1.526434 

FIRST 

HIT 

REAL" 

1.300000 

XMAG* 

.850000 

zeta* 

.836970 

OMEOA" 

1.553222 

LAST 

REAL" 

REAL" 

l.AOOOOO 

l.AOOOOO 

IHAO« 

IHAO« 

•050000 

•700000 

ZETA* 

zeta* 

•999363 

•894427 

OMEOA" 

OMEOA" 

1.400893 

1.565248 

FIRST 

HIT 

REAL" 

l.AOOOOO 

XHAG« 

.750000 

ZETA* 

•881480 

OMEOA" 

1.580238 

HIT 

REAL" 

1.400000 

IHAG« 

•800000 

ZETA* 

•868243 

OMEOA" 

1.612452 

LAST 

REAL" 

'REAL" 

1.500000 

1*500000 

IMAO* 

VlMAGa 

.050000 

•600000 

zeta* 

ZETA* 

•999445 

.928477 

OMEOA" 

OMEOA" 

1.500833 
1.615569 

FIRST 

HIT 

REAL" 

1.500000 

xnao* 

•650000 

zeta* 

.917556 

OMEOA" 

1.634778 

HIT 

REAL" 

1.500000 

xnag* 

*700000 

zeta* 

•906183 

OMEOA" 

1.655295 

LAST 

REAL" 

REAL" 

1.600000 

1.600000 

xmag« 

XHAQ« 

.050000 

.450000 

zeta* 

ZETA* 

.999512 

.962651 

OMEOA" 

OMMA" 

1.600781 

1.662077 

FIRST 

HIT 

REAL" 

1.600000 

inag- 

•SOOOOO 

ZETA* 

•954480 

OMEOA" 

1.676305 

HIT 

REAL" 

1.600000 

ZHAO» 

.550000 

zeta* 

.945687 

OMEOA" 

^691892 

LAST : 

REAL* 

REAL" 

1.700000 

1.700000 

IMAQa 

imag« 

•050000 

.OSOOOO 

zeta* 

ZETA* 

•999568 

•999568 

OMEOA" 

OMEOA" 

1.700735 

1.700735 

FIRST 

HIT 

REAL" 

1.700000 

XHAG« 

•lOOoOO 

ZETA* 

•998274 

OMEOA" 

1.702939 

HIT 

i REAL" 

r<l700000 

IMAO« 

•150000 

zeta* 

•996130 

OMEOA" 

1.706605 

HIT 

REAL" 

"1*700000 

IMAg> 

*200000 

ZETA* 

•993151 

OMEOA" 

1.7U724 

HIT 

REAL" 

1.700000 

XHAG- 

•250000 

ZETA* 

•989359 

OMEOA" 

1.718284 ' 

HIT 

REAL" 

"1.700000 

xhag« 

.300600 

ZETA" 

.984766 

OMEOA" 

1.726288 

HIT 

REAL" 

'1.700000 

IMAG* 

•350000 

ZETA* 

•979457 

OMEOA" 

1.735655 

LAST 

POLEi 

oviR'ia • 










. CJ' 

"V” m<L-yL 


*1000 

• 0714 

•>.2285 

*2000 

• 2812 

•.6796 

.3000 

• 6136 

• 1.8648 

.4000 

U 0314 

• 3.0752 

• 5000 

1*4981 

• 4.4106 

.6000 ' 

1*7984 

- 5.7941 

.7000 . 

1.7780 

- 7.1734 

.6000 

1*9690 

- 8.5463 

.9000 

1*9261 

• 9.9803 

1 . 0000 ' 

1*5797 

- 11.3373 

I.IOOO 

*9278 

• 12.6121 

1.2000 

*0358 

• 13.8057 

1*3000 

- 1.0024 

• 14.9219 

[ 1.4000 

« 2*1026 

- 1 S . 96 S 7 

11.9000 

- 3.0691 

- 16.9426 

1 1.6000 

- 3*6366 

• 17.8581 

1.7000 

- 4,2167 

• 18.7174 

1*8000 

- 4.7988 

• 19.5252 

1.9000 ' 

- 5*2337 

• 20.2860 

2.0000 

- 5*6656 

- 21.0038 

2.1000 

- 6.0762 

• 21.6820 

2.2000 

• 6*4376 

• 22.3241 

2.3090 

- 6*7923 

- 22.9328 

2.4000 

- 7*1400 

- 23 . 5 il 0 

2.9009 

- 7.4805 

• 24.0609 

2.6000 

- 7*8136 

- 24.5848 

2.7008 . 

- 8*1393 

• 25.0846 

2.6000 • 

- 8*4575 

- 25.5621 

2.9000 ' 

- 6*7684 

• 26*0189 

3.0000 

- 9*0720 

- 26.4565 

, 3.1000 

- 9.3684 

• 26.8762 

! 3. 2000 

- 9*6579 

- 27.2793 

[ 3.3000 

- 9*9405 

• 27.6668 

3.4000 

- 10.2165 

- 28.0398 

3.8000 

- 10*4863 

• 28.3993 

3.6000 

- 10*7495 

- 28.7460 

1 3*7000 

- U .0069 

- 29.0809 

3.8000 

- 11*2585 

- 29.4045 

3*9000 

- a . 5044 

- 29.7177 

4.0000 

- 11*7450 

- 30.0210 

4.1000 

- 11*9803 

- 30.3150 

4.2000 

- 12.2106 

- 30.6002 

4.3000 

- 12.4361 

• 30.8771 

4.4000 

- 12*6969 

• 31.1462 

4*5000 

- 12*8732 

- 31.4078 


\ 


e 

• 2999 

-7.1978 

-20.9959 

1.1608 

-16.5785 

-41.0431 

2.4784 

-22*3371 

-60.0069 

4.1067 

-30*6910 

-77.8711 

5.8688 

-39,8816 

-93,7696 

7.5525 

•50.1576 

•107,5548 

8.9514 

•58.8466 

-119.3672 

10.5153 

-66^5168 

' -129.4850 

11.9064 

-73*9336 

-138,2087 

12.9170 

•81*0828 

•145.8059 

13.5400 

•87*9580 

-152.6957 

13.5415 

•96*5560 

-158.4509 

13.9195 

-100*8T75 

•163.8056 

13.8631 

•106*9263 

•168.6641 

13.8735 

-112*7096 

-173.1074 

14.2215 

-118*2316 

-177.TB54 

14.5007 

-123*5052 

•182.2910 

14.7265 

•128*5392 

•186.3371 

15.0523 

•133*3463 

•190.0116 

15.3382 

-137*9316 

•193.3801 

15.6058 

-162*3115 

•196*4926 

15.8865 

-146*6953 

-199.3871 

16.1405 

-150*4933 

-202.0949 

16.3710 

•154;3l59 

•204.6387 

16.5804 

-157*9T26 

-207.0391 

16.7712 

-161*4729 

-209.3110 

16.9653 

•164.8257 

-211.4700 

17.1046 

-168^0392 

•213.5250 

17.2505 

•171.1214 

•215.4881 

17.3845 

-174.0T97 

•217*3616 

17.5078 

-176*9211 

•219.1585 

17.6214 

-179.6523 

-220.8829 

17.7263 

-182*2792 

•222.5602 

17.8233 

-184^8077 

-224.1353 

17.9131 

-197*2*32 

-225.6724 

17.9965 

-189*5908 

•227.1553 

18.0740 

-191.8552 

•220.5875 

18.1461 

•194.0407 

-229.9722 

18.2133 

-196*1517 

-231*3121 

18.2760 

-198*1919 

•232.6100 

18.3347 

-200*1650 

-233.8681 

18.3896 

-202*0744 

-235.0887 

18.4610 

-203.9233 . 

-236.2738 

18.6893 

-205*7148 

-237.4253 

18.5366 

-207.4516 

•238.5448 


IS. mri 

4!««4646i 
S7.66«T; 
47.1601 
, 93.8678 
57.3972 
60.6228 
62;8682 
64.2792 
64.7230 
64.9377 
63.8949 
62.9281 
61*7378 
60.3990 
' 59.9938 
98 * 7898 :, 
67.79791 
56 . 6673 ' 
55*4487 
94 . 1809 , 
92 . 8919 ; 
91.6012 
50.32281 
49.0669 
47.6388 
46.6444 
49.4858 
44.3648 
43.2819 
42.2374 
41.2306 
40.2610 
39*3276 
38 . 4291 ; 
37.5645 
36.7323 
35.9314 
35.1605 
34.4181 
33 . 7031 ; 
33.0143 
32.3505 
31 . 7105 ; 
31*0932 



6.6000 

•13.0853 

•31.6625 

18.5773 

-209.1365 

•23986341 

6.7000 

•13.2931 

-31*9105 

18.6176 

-210.7718 

•240t694S 

6.8000 

•13.6970 

-32.1523 

18.6553 

-212*3000 

•241. 7275 

6.9000 

•13.6970 

-32.3880 

18.6910 

-213.9031 

-242.7344 

s.oooo 

-13.8933 

-32.6181 

18.7267 

-215*6036 

-243.7164 

S. 1000 

•16*0861 

-32.8627 

18.7566 

-216.8628 

-244.6746 

5.2000 

-16.2756 

-33.0622 

18.7868 

•218*2830 

-245.6102 

5.3000 

-16,6616 

-33.2768 

18*8156 

-219*6659 

-246.5240 

5.6000 

•16.6662 

-33.4867 

16.8625 

-221.0130 

-247.4171 

5.5000 

•16.6239 

-33.6921 

18.6682 

-222.3260 

-246.2903 

5.6000 

-15.0006 

-33.8933 

18.8927 

-223.6062 

-249.1445 

5.7000 

-15.1765 

-34.0906 

18.9159 

-226.8550 

-249.9804 

5.8000 

•*lS,34Sb 

-34.2836 

18.9380 

-226.0736 

-250.7988 

5.9000 

-15.5160 

-34.4730 

18.9590 

-227.2638 

•25U6004 

6.0000 

-15.6799 

-34.6589 

18,9791 

•228.6261 

-252.3859 

“BTTOffO 

-15, 8632' 

■ -34,8614“ 

18.9992 " 

-229.3619 ■ 

-253.1558“ 

6.2000 

-16.0062 

-35.0206 

19.0166 

-230.6721 

-253*9108 

6^3000 

-16.1628 

-35.1966 

19.0338 

•231.7579 

-254.6514 

6.6000 

-16.3192 

-35.3696 

19.0506 

-232.8200 

-255.3782 

6.5000 

-16.6736 

-35.5397 

19.0663 

•233.8595 

-256>0917 

6.6000 

-16.6256 

-35.7070 

19.0816 

-236.8771 

•256*7923 

6.7000 

-16,7755 

-35.8716 

19.0961 

-235*8738 

-257*4805 

6.8000 

-16.9235 

-36.0336 

19.1101 

-236*6501 

-258*1567 

9*9000 

-17.0697 

-36.1931 

19*1236 

-237.8070 

-258*6214 

7.0000 

-17,2160 

-36.3502 

19,1363 

•238.7650 

-259*4749 

7.1000 

• -17.3566 

•36.5050 

) 9* I486 

-239*6668 

-280*1177 

7.2000 

-17.6972 

-36.6575 

19.1606 

-260.5671 

-260.7500 

7*3000 

-17.6362 

-36*5079 

19.1717 

-261*6526 

-261*3721 

7.6000 

-17,7736 

•36.9562 

19*1826 

•262*3216 

-261*9845 

7*5000 

-17.9096 

-37.1025 

19*1931 

-263*1765 

-262*5875 

7.6000 

-18.0637 

-37.2668 

19.2032 

-266*0i23 

-263*1813 

7.7000 

-18.1766 

•37.3893 

19*2129 

-266.8353 

-263*7661 

7.8000 

-18,3077 

-37.5299 

19.2222 

•265;6460 

-264.3424 

7.9000 

-18.6375 

-37,6687 

19*2312 

-266*6388 

-264*9103 

8.0000 

-18.5660 

-37.8059 

19*2399 

-267*2202 

-265.4700 

8.1000 

-18.6931 

•37.9613 

19*2682 

-267*9888 

-266*0219 

8.2000 

-18.8189 

-38.0752 

19.2563 

-268*7662 

-266.5661 

6*3000 

-18,9636 

-38.2075 

19.2661 

-269»6i76 

-267*1026 

8.6000 

-19.0667 

-38.3383 

19.2716 

-250*2192 

-267*6323 

8.5000 

M -19.1888 

•38.6676 

19.2788 

-250*9392 

-268*1548 

8.6000 

-19,3097 

-38.5955 

19.2856 

-251*6680 

-266*6704 

8.7000 

T -19.6296 

-38.7221 

19.2926 

-252*3659 

-269.1793 , 

8.8000 

-19.5681 

-38.8672 

19.2992 

-253*0332 

-269*6817 

8.9000 

-19,6656 

-38.9711 

19.3055 

-253,7103 

-270.1778 


30**976. 

29,9227 

29,3676 

2B.B313 

2B.3130 

27*8nB 

27.3271 

26*BBB1 

26.6060j 

23.9663 
25.53B3 
2B.1256i 
26.7230 
26.3366 
23. 9898 
“2575939 
23.2387 ' 
22.8936 
22.5582 I 
22.2322. i 
21.9152 
21.6067 
21.3066 
21*0165 
20.7300 
29*6529 
20.1829 
' 19.9197 
19.6632 

19.6130 I 
19.1689 
18*9308 i 
18.6986 
18.6716 
18.2698 
18.0336 
17.8219 
17.6152 

17.6131 
17*2156 
17.0226 
16.8336 
16.6685 
16.6675 



